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We introduce a class of quantum channels with correlations acting on pairs of qubits, where the 
correlation takes the form of a shift operator onto a maximally entangled state. We optimise the 
output purity and show that below a certain threshold the optimum is achieved by partially entangled 
states whose degree of entanglement increases monotonically with the correlation parameter. Above 
this threshold, the optimum is achieved by the maximally entangled state characterizing the shift. 
Although, a full analysis can only be done for the 2-norm, both numerical and heuristic arguments 
indicate that this behavior and the optimal inputs are independent of p > 1 when the optimal 
output purity is measured using the p-norm. 
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I. INTRODUCTION 

In usual memoryless channel model successive uses of 
the communication line have the same noise [l[ and can 
be described as a simple tensor product of channels. Re- 
cently, there has been some interest in studying the be- 
havior of these channels with correlations 0)03] since such 
channels might be regarded as a small first step in study- 
ing the much more complex issue of channels with mem- 
ory @, H, [1, 0, H, [§]. They can also describe multiple 
access channels. Thus, one is led to consider a scenario 
in which with probability 1 — fx qubits encounter only 
uncorrelated tensor product noise, but with probability 
fi they experience correlated noise. This situation can be 
modelled by a channel of the form 



$ = (1 - /i)* ® * + fiT c 



(1) 



where all the correlations are in the map r corr , which 
need not be a channel itself. Several papers 0, y] have 
studied specific examples of this type for which there is 
a critical value /i c below which the optimal input is a 
product state and above which the optimal input is max- 
imally entangled. Recently, Daems [10| showed that this 
is always the case for a class of channels called "product 
Pauli". This is not, however, the most general channel 
of the form (TT|), even for a pair of qubits. 

In this note, we consider a class of channels with corre- 
lations which exhibit quite different behavior. The map 
we use for r corr in |T]) is extremely simple; it is simply 
the non-unital channel which maps every input to a fixed 
maximally entangled state. We study the optimal out- 
put purity, rather than the channel capacity. For our 
channels, the entanglement of the output input increases 
continuously with /j, until it reaches a critical /z c , after 
which the optimal input is always achieved with a maxi- 
mally entangled state. 

Although one is ultimately interested in the effect of 
correlations on various types of channel capacity, we con- 
sider here only the optimal output purity. For channels 
with some covariance properties [ll| , one can make an ex- 
plicit connection between the classical capacity [l2| and 



the optimal output purity as measured by the minimal 
output von Neumann entropy. However, that need not 
hold in general and is not true for the channels consid- 
ered here. Moreover, the conjecture in [l3| about capac- 
ity achieved with maximally entangled states depends on 
the precise form of the channels and is not relevant here. 

Roughly speaking, one expects inputs whose outputs 
are close to pure states to be the least corrupted; how- 
ever, one can have very noisy channels which map all 
inputs to a region close to a fixed pure state with lit- 
tle correlation with the input. Nevertheless, the optimal 
output purity as measured by either the minimal output 
entropy [TJ, [l5[ or the maximal output p-norm [l6[ is 
of some interest and has been studied extensively. This 
is, in part, due to the Shor equivalence 17[ established 
between the conjectured additivity of minimal output 
entropy and other long-standing additivity conjectures. 
The hope that the additivity of minimal output entropy 
could be proved by showing that the maximal output 
p-norm is multiplicative, at least near ji « 1, as conjec- 
tured in [l6j , was recently shattered by the discovery jl8| | 
of counterexamples to the latter for all p > 1. Neverthe- 
less, the additivity conjectures remain open and even the 
multiplicativity conjecture is known to hold in the re gion 
1 < p < 2 for certain classes of channels [l9l . I20I l2ll J22| 
and in other situations for p = 2 [23|, [24|. Thus, the op- 
timal output purity remains an object of some interest. 

For p > 1, the p-norm of a state 7 is given by the 
expression 



l7llp = [Tr(7 P )] 



(2) 



One sometimes uses instead, the Renyi entropy [25[ 

^(7) = ^ In Tr(<f) = 1 ^]n\\< r \\ p , (3) 

which is known to converge to the usual von Neumann 
entropy as p — > 1. The maximal output p-norm of a CPT 
map $ denoted is the supremum of ||3>(7)||p over 

all input density matrices 7, i.e., 



=sup ||$(7)||p 

7 



(4) 
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The minimal output entropy and Renyi entropy are sim- 
ilarly defined as 

S p ($)=MS p [m)], (5) 

7 

and it is natural to refer to states which achieve the op- 
timum in ^ or ([5]) as optimal inputs. 

This paper is organized as follows. In Section [IT] we 
describe the class of channels we study and show how 
covariance properties can be used to reformulate the op- 
timization problem. In Section IIIII we use unitary trans- 
formations to simplify the problem, solve it exactly when 
p = 2, reduce the general case to the analysis of a single 
parameter, and report numerical work which supports 
the conclusion that the optimal inputs are independent 
of p. In Section IIV[ which can be skipped on first read- 
ing, we analyze the behavior of the output eigenvalues 
under certain small perturbations. Although this analy- 
sis does not yield a proof for p ^ 2, it does support our 
conjectures and yield a proof for p — oo. In Section fVlwc 
summarize our conclusions in the form of both conjec- 
tures and theorems, and summarize the evidence for the 
former. These can be reformulated as statements about 
the trumping relation [H, HI HI HJ [!j which plays an 
important role in entanglement catalysis. 

In addition to the isomorphism C4 ~ C2 £3>C2, there is 
also an isomorphism between vectors in C4 and matrices 
in M2, and a straightforward way to make this corre- 
spondence using the Pauli matrices. We describe this in 
Appendix [AJ Although the results are straightforward 
and/or well-known, they play an important rule and it is 
useful to describe them in a fixed notation. 



where 

|0) = (tft ® yt) |ft> = (I ® yt u) |/3 ) (10) 

and the second equality used (|A3[) . Since unitary trans- 
formations do not affect eigenvalues, 

® V ) R ^ ® ^ f )] Hp = \\^^AR)\\v (11) 

when |/3) is given by (I10|) . Because the p-norm is convex, 
it suffices to consider the optimization j3]) for pure states 

We now define an equivalence relation on pure states 
W) e C 4 by 

l^i) = \*h) <* 3 unitary U, V : U <g> V\i/>i) = |^ 2 ) . 

(12) 

Since all members of a given equivalence class are re- 
lated by local unitaries, they have the same entangle- 
ment. In the matrix picture described in Appendix [Al 
each equivalence class is characterized by its singular val- 
ues. Moreover, we can characterize each equivalence class 
by its entanglement as measured by the entropy of its 
reduced density matrix (or by replacing h in (|A8j) by 
another function strictly monotone on [0, 1]). In particu- 
lar, we find it useful to use the so-called "linear entropy" 
E = 2(1 — Tr7 2 ) in numerical work. 

We now let l^) denote an equivalence class or, more 
properly, a representative of the class with properties to 
be specified in the next section. Then it follows from l(TTj) 
that 



II. CHANNEL DEFINITION AND PROPERTIES 

Let denote the qubit depolarizing channel [12] 

* A ( 7 ) = (l-A)±J(Tr 7 ) + A 7 (6) 

with A G [— 1 , 1] . From this we generate a correlated 
two-qubit channel of the form (JTJ) whose action is 

*/s,p,A(i2) - (1 - A*)(*A ® + ^ (Tri?) |/3) (0| ,(7) 

where < p < 1 and \j3) is a fixed maximally entangled 
state of the two qubits. Although is well-defined for 
any matrix in M4, we are interested in the case of density 
matrices, for which R > and Tr R = 1. We will exploit 
the covariance property of the depolarizing channel, i.e. 

(8) 



which holds for any unitary U matrix in Mi- We find a 
relationship between the channel obtained using the Bell 
state |/? ) = -75(100) + |11)) and any other maximally 
entangled state by observing 

<S>p ^x[(U®V)R(U^ ®V^)] 



sup^ II^o.m.aCIV'X^DIIp 

= supsup[|^ 0)M ,A[(t/®y)|^)W(C/ t <s> V^)] h„ 



SUpSUpH^^AdW^DIIp 

£ 



(13) 



with |/3) maximally entangled. 

The observations above allow us to draw several con- 
clusions 

a) First, ((lip implies that \ip) is an optimal input for 
&p,ti,\ if an d only if U (g> V'lV') is an optimal input 
for $/3 0l/tl A- Therefore, 



^p(^/3,m,a) = ^(^ ,m,a) 



(14) 



The same conclusion can be reached by reversing 
the roles of |/3o) and |/3) in (|T5)| . Thus, it is sufficient 
to study ®0 o<fl< \. 

b) If the optimal |?/>) in the reformulation (|13[) is 
unique, then the set of optimal inputs for §p 0ifl ,x 
is a subset of {|^>) = (ifi ® V rt )|V') : ?7, y unitary}. 
Thus, we expect that, excluding some trivial cases 
(A = or p, = 1), all optimal inputs for a given 
(9) channel have the same entanglement. 



yt) 
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c) The singlet state \02) = f2 1 /3o } satisfies the covari- 
ance condition U <8> U\(3 2 ) = I ^2) for all unitary 
U G M2. Therefore, if we modify the equivalence 
relation (TTZ|) by restricting to U = V, 



fj,,\) 



sup 



I^/3 2 ,m,a(IV')(V'I 



(15) 



where \ip) denotes a representative of the modified 
equivalence class. 

Before performing the optimization, we mention some 
trivial cases, always using the channel with \(3q) for the 
shift. When A = or [i = 1, all inputs have the same 
output, namely \{\ — + /j,\f3o)(Po\- We will see that 
these are the only cases in which inputs with different 
entanglement are optimal. When fj, = 1, the optimal 
input \Po) is unique and the optimal output is also |/3o)- 
Henceforth, we will assume that A € (0, 1) and fi G [0, 1). 
When /i = every product state is an optimal input; at 
the end of Section flllDI we prove that, as expected, all 
other states are non-optimal. 



III. OPTIMIZATION 

A. Simplifying the input 

After taking into account the normalization condition 
and irrelevance of an overall phase factor, the optimiza- 
tion problem for v p ($/3 0tl j, t \) involves three complex, or 
six real, variables. The covariance used to obtain (fl3|) 
reduces this to four real variables, one for ip and three 
for Moreover, the additional symmetry noted in (c) 
above allows an immediate reduction to three real vari- 
ables. We use a different approach; the reduction from 
four to three real variables is obtained following (J22). In 



Section UlIDl we find that it would suffice to analyze the 
dependence on a single variable with the others fixed. 

Using the notation of Appendix [5] for the maximally 
entangled Bell states, we can write an arbitrary state 
tp E C 4 as 

\i>) =^> fc |/? fc ) =J2a k (I®a k )\/3 ) = (I (g> A)\Pq) , (16) 



where A — J2k ak<Tk e ^2- Moreover, we can use 
the SVD to choose unitary matrices U, V so that (U ® 
V)\ijj) = (I <g> VAU T )\Po) with VAU T diagonal and pos- 
itive, as discussed in Appendix El It will be convenient 
to write the corresponding state (|A5|) using an angular 
variable 9 £ [0, so that 

(U®V)\il>) = \il> e ) = 75(cosf + sin §)!#,) 

+ -^(cos|- S m|)|/3 3 ) , (17) 

which implies 

\^e)(ipe\ = \ I® I + <Jz + cos 9 (a z <g> I + I ® a z ) 



+ sin6(a x ® a x — a y ® tXj 
B. Computing ||$ /3 , M ,A(|V'e)(V'e|)llp 



(18) 



Since ^(ofc) = Ao^ it is straightforward to see that 
(TTSJ) implies 



(*A®*A)(|^9><^fl|) 



(19) 



= J I®J+ \ 2 <r z (g)(T z + cos 9\{a z ® I + I ® a z ) 
+ sin 9 A 2 ia x ® a x — o y <g> a y ) , 

which can be written as 



(*A®#a)(|^)M = \ 



/1 + A 2 + 2A< 





Now writing V^U = ( a ^ !L ) will) <i 



V 2A 2 sin6> 
2 + \b\ 2 = 1, we find 





(1-A 2 ) 





2A 2 sin6> 

(1-A 2 ) 

l + A 2 -2Acos( 



(20) 



= (/3 |(FtC/)t = ( / 9 |(J a b 



= (a —b b a) 



Using (f20|) and (j2Tj) we find that &f3, l _ l ,\(\i>e)(ipe\) is given by 

T + A 2 + 2Acos6> 2A 2 sin6> 

1 -// I 1-A 2 

1-A 2 

2A 2 sin(9 l + A 2 -2Acos( 



2 



/|a| 2 — ab ab a 2 \ 

1 -ab \b\ 2 -b 2 -ab 1 

ab -t \b\ 2 ab 

\ a 2 — ab ab |a| 2 / 



(21) 
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It is evident that (0 b b 0) T is an eigenvector of both matrices above. This suggests that we make a simplification 
exploiting the block structure in the left matrix using a basis which includes the known eigenvector. We act first on 
both matrices with the permutation matrix P which exchanges the 2nd and 4th rows and columns, and then make a 
unitary transformation which preserves the block structure and achieves a partial diagonalization. Thus, we replace 
each matrix above by W^P( )PW where 



P 



(1 0\ 

1 

10 

Vo i o 0/ 



f l fa a 
M la —a 



\ 







1 

[61 





b 
b 



(22) 



Although these unitary transformations do not preserve entanglement, they do not affect the eigenvalues of the output. 
After introducing the shorthand c v = cos^, s v = smip,S — 2Asin#,C = 2Acos#, and M M = 4/i/(l — /«), we find 
that the transformed output density matrix is 



/l + A 2 



C 



a AS* 



V M M VN 2 (1-H 2 ) 



C + is<f 
1 + A 2 - 





xs 



, xs 






1-A 2 



1-A 2 






-M M (1 



0/ 



(23) 



We have now reduced the optimization to a problem in 
three real variables, sin#, ip, \a\. 



C. Optimization for p — 2 



It follows from (|23[) and some elementary algebra that 



(24) 



= C-^Y{mI + 2M42A 2 |a| 2 + (1-A 2 ) + c v XS] 
+ 4(1 + A 2 ) 2 -2S 2 (1-A 2 )} 

Since only the coefficient of Af M includes any dependence 
on a and M M > 0, it follows that the 2-norm is largest 
when this coefficient is largest (see Fig. [T]). This occurs 
when \a\ — 1, and c v = 1. After making these choices we 
find 

|| || 2 = (i^) 2 ^/ 2 + 2M AI (l + A 2 + A5) 
+ 4(1 + A 4 ) + 2(C 2 + A 2 ^ 2 )] 
= ( i 4 £ ) 2 [M 2 + 2A/ AI (l + A 2 + 2A 2 sin6») 
+ 4(1 + A 2 ) 2 -4A 2 (1-A 2 )sin 2 6»] (25) 

which can be regarded as a quadratic function of sin# 
whose optimization is straightforward. Since we obtain 
an equivalent problem whenever the optimum is achieved 
with \a\ = 1, we give the details in the next section. 

D. Consequences of |o| = 1 optimal 

We now describe the conclusions one can reach if the 
optimal output p-norm is achieved when \a\ = 1. We will 



show that when this is the case, the optimal inputs are 
the same as for p — 2. 

When \a\ = 1, the output density matrix (j23f is block 
diagonal. Two of its eigenvalues are |(1 — /t)(l — A 2 ) and 
the remaining two are 



i(l- M )(l-A 2 ) + i/i 



(26) 



± V)\[\M~1 + 4A2 - (! - A2 ) 52 + M^XS . 

The output p-norm will be largest, and output entropy 
smallest, when the quantity under the square root is 
largest. Therefore, one should choose c„ = 1. This im- 
plies that V^U — I or, equivalently, that V = U. 

The optimization problems for both |25|) and (|26|) then 
reduce to maximizing a function of the form 



4(1 - A 2 )A 2 sin 2 i 



2A^M AI sin<9 



constant (27) 



which is quadratic in sin# > when A £ (0, 1). This is 
largest when sin 6* = min{l, }■ We distinguish 

two situations characterized by the threshold value \i c = 

1 —A 2 

' corresponding to the boundary sin = 1. 



2-A 2 



• Below Threshold: If /t < /t c then the maximum 



value of (1271) is achieved for 4 sin 9(1 — X z 



fi(l - /t c 



M„ 



4^^, or, equivalently 



l—/< 1-A 2 



(28) 



Moreover, since V = U , the optimum is achieved 
with a family of input states of the form 

IVm) = [V T ®V^)\^ e ) , (29) 

with 9 as in (|28p. \ipg) as in ([T7]), and V an arbitrary 
unitary. Using (|A7|) . one finds that the reduced 




FIG. 1: Plot of the 2-norm for the channel $^ 0iMi a- The 
continuous curve represents the threshold boundary. For n = 
1 or A = 1 the norm is 1: in the former case the channel sends 
every input state into |/?o)(/3o|; in the latter, it transforms 
the input (i/>| into a mixture of \ip){ip\ and |/3o)(/3o| so that 
choosing = \Po) gives a pure output. 
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FIG. 2: Entanglement of the optimal input states which max- 
imize the 2-norm for the map &p 0tlJ ,,\. For \x € [0,/j. c ) the 
entanglement of the optimal inputs is a strictly monotone 
function of fx; for /i > /i c the optimal input is the maximally 
entangled state |/3o) - 



density matrix of \ipg) is 

76 = \[I + 008 61(73] (30) 

and the entanglement of any optimal input is sim- 
ply h(cos9) where h is the binary entropy (|A8j) . or 
E = sin 2 9 when the linear entropy is used (see Fig. 

• At or above Threshold: For fj, > fi c the maxi- 
mum value of (|27[) is achieved for sin 9 — 1. Then 
(fT7| gives \ipe) = |/?o) and the diagonal matrix Dg 
is simply the identity so that DgV = V'V = I. 
Thus, there is a single optimal input, the maximally 
entangled state |/3o)- 

Note that for /i = the value of \a\ is irrelevant and 
the optimization gives a special case of (1271) with the op- 
timum achieved for 9 — 0, consistent with (f28|) . Because 
there is no dependence on a in this case, this result holds 
for all p and for all \f3) in (fT3|) . Hence, V'U is arbitrary 
and we recover the expected result that all product in- 
puts are optimal. Because 9 ^ can never yield a state in 
the equivalence class with no entanglement, we also find 
that a state which is not a product can not be optimal. 

E. Numerical results 

It is not easy to perform an exact analytical analy- 
sis of the output p-norms for p ^ 2. However extensive 
numerical studies of optimization were carried out using 
the equivalent Renyi entropy [IH, for over 2000 pairs of 
randomly chosen value of \i and A. In all cases, we found 
that the input states which are optimal for p — 2 are 



also optimal for p > 1. In Fig. [3l we show typical nu- 
merical results of our findings by comparing the minimal 
Renyi output entropies ([5|) for randomly chosen inputs 
with that of optimal inputs for p = 2, which lie on the 
bottom curve. In all cases the output Renyi entropy was 
larger than the expected minimum. 



IV. PARTIAL ANALYSIS WITH HEURISTICS 

A. Eigenvalue behavior from characteristic 
polynomial 

In this section, we attempt to show that the optimal 
p-norm is attained when |et| = 1 be showing that it in- 
creases monotonically with \a\. Although our argument is 
incomplete and can not exclude fluctuations under some 
conditions, it does provide additional support for this 
conjecture. To analyze the behavior of the eigenvalues of 
(|23[) as a function of the parameters \a\ and if, we use the 
characteristic polynomial of the output density matrix to 
estimate the effect of small changes after reduction to a 
3-dimensional problem. 

Since ^(1 — p)(l — A 2 ) is clearly an eigenvalue of (|23j) . 
we are left with the eigenvalue problem for the 3x3 
matrix 

A = (l + A 2 )/+ (31) 

/ c v XS + Mt,\a\ 2 C + i Slfi \S M^\a\ 2 {l- \a\ 2 )\ 

C - is v AS -c v \S 

\A/ M V>| 2 (1 " M 2 ) -2A 2 + M(1- \a\ 2 )) 

where C, S, <p and are defined as in (|2"3l) . 
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FIG. 3: Scatter plots of the output Renyi entropy 
S'p[ < 1 , /3 ,h,a(|^){V , I)] as a function of p for randomly chosen in- 
puts compared to that for the conjectured optimal input. The 
horizontal line corresponds to the maximum possible Renyi 
entropy of In 4. Top: A = |, [i = | < /i c 



old. Bottom: A = i, /i 



| below thresh- 



-pf above threshold. 



The characteristic polynomial for A — (1 + A 2 )/ is 

%) = ~( 3 + (M^-2\ 2 )C 2 +[C 2 + \ 2 S 2 + M^\a\ 2 (2\ 2 
+ c tp XS)](+(C 2 + X 2 S 2 )(2X 2 - M) 
+ M>| 2 (C 2 + \ 2 S 2 + 2\ 3 c ip S) 



2 (2X 2 

M t ,\a\ 2 (C 2 + X 2 S 2 + 2X 3 c ip S) . 



i?o(C) + M>| 2 (2X 2 + c v XS)( 



(32) 



This defines R (()- Tncn with a; = C + 1 + A 2 , the char- 
acteristic polynomial for A is 



R(x) = R (x - 1 - A 2 ) + M^\a\ 2 (2X 2 + Cip XS)x 



(33) 



+ (l-A 2 )(2A 2 -c v AS*-5 2 ) 
where we used C 2 + S 2 = AX 2 . 

B. Optimization for 1 < p < 2 



We first assume that a, S are fixed and consider the 
optimal value of c v ~ cos tp by examining the effect of 



the change c ¥ 



P(x) = R(x) given by (|3"3"j). and Si 



e . We apply Lemma [6] with 



eMn\a\ 2 XS and 



S2 = — e M p |a| 2 AS'(l — A 2 ) to find that the eigenvalues of 
A 



1 



[«fc - (1 - A 2 



(34) 



with fc, to, rt distinct. It then follows from Lemmas [7] and 
[5] using the expressions (|B8|) and (|B10|) that 



| V ||P_> llvF + eAfJa| 2 A5 



(.-Dd-A^l^-^ J i +0(e 2 )(35) 



where «i > «i > 1*2 > «3 > f3 and similarly for iik . Then 
the quantity in square brackets [ ] > is positive for 
1 < p < 2. Since AS* = 2A 2 sin 6* > (by our assumption 
9 e [0, -|]), \\A\\p increases as c v goes from —1 to +1 and 
is thus largest when c v = 1. 

To study the effect of changing \a\, first observe that 



when 



1, 2A 2 



c v XS 



S 2 = 2A 2 (1 + sin0)(l 



2 sin 6*). If we insert this in (f3"3")> . we can apply Lemma [5] 
with P{x) = R{x), Si = e2A 2 A/ M (l + sm6) and S 2 = 
e2X 2 M fl (l + sin0)(l - A 2 )(l - 2sin6») to conclude that 



v k + e 2A 2 M M (1 + sin(9) 



1 



\Vk 



{vk - v m )(v k - v n ) 
+ (l-A 2 )(l-2sin6»)] (36) 

We again apply Lemmas [7] and [8] using the expressions 
(|B8)l and (|BT0)) to conclude that 



\\ p + eB 



p(vl^-vl^) + (p~l){l-X 2 ) 



(l-2sm6){v p 1 - 2 ~vP- 2 )] +0(e 2 ) (37) 

with Vk,Vk constrained as above and B = 2A 2 M M (1 + 
sin 8) - g- > 0. However, v p ~ 2 < v p ~ 2 because p — 2 < 
0. Therefore, we can only conclude that the quantity in 
brackets is positive when sin6> > i. Otherwise, (|3"5|) has 
two competing positive and negative terms. 

When (i > ^Z\'i 1 the optimum in (|28[) satisfies sin 9 > 
■i. Even if the p-norm does not increase monotonically 
with \a\ for small values of 9, it seems likely that the 
optimum is still achieved when \a\ = 1. 



C. Optimization for p > 2 



When p > 2, the expression (|35|) contains competing 
terms and we cannot reach a definite conclusion about 
the effect of c M 1— > c„ + e . However, both (l34l) and (l36l) 



still imply that the largest eigenvalue increases under the 
changes c v 1— > + e and |a| 2 1— > |a| 2 + e with e 
Moreover, if one fixes <p and considers the change 
\a\ 2 + (Six + S 2 ) with 



> 0. 

,|2 ^ 



Si = e2M M A 2 (l 



(38) 



e2M M A 2 (l - A 2 )(l - cosine - 2 sin 2 9) (39) 
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then Lemmas [6] and [7] imply that 

||A||£ » WAWV + eBpil + CvSmeM- 1 -**- 1 ) 
+ eB(p-l)(l - A 2 )(l - c v sin0 

- 2sin 2 0)(i;r 2 -«r 2 ) + O ( e2 ) (4°) 

where = 2X 2 M^p. When sin^ < | this implies that 
|| A ||p increases with |a| 2 . This is sufficient to show that 
when /i < g at least a local optimum is achieved when 
c v = \ and is given by (J25J). 

D. Additional heuristics 

One can apply the MVT again to (|57|) using t>2 , i>2 to 
denote the mean values. Under the assumptions Vk ~ Vk 
and i>2 > (1 — A 2 ) the term in square brackets [ ] in (I57| 
is 

« (p - 1)^(1 -A 2 )2[Cp-I)(l- sin 0)+sin#] > 

(41) 

when 1 < p < 2. For p > 2, a similar analysis beginning 
from (|40p gives a change in ||A||*j approximately propor- 
tional to 

(p-2)(f-sin 2 6») + (f + C[p sm6) > 0. (42) 

This approach can even be applied to the entropy to show 
that when |a| increases the output entropy decreases un- 
der the assumptions above. We omit the details. 

Since v\ > i\ > > V3 > V3, one knows that ii 
is larger than the smallest eigenvalue of A and proba- 
bly close to the second largest, which one expects to be 
> 1 — A 2 . Thus, the assumptions above are reasonable. 
However, this is far from the desired proof that the out- 
put p-norm increases with \a\ and is thus optimal for 
\a\ = 1. For p > 2, these heuristics are more convinc- 
ing because errors from these approximations are better 
controlled . 

V. CONCLUSIONS 

A. Main results 

Conjecture 1 Let $p ^ x be a channel on M4 as defined 
in (|7|) with fi, A G (0, 1) and let p c — \Z_\2 ■ Then 

i) For < p < p c the maximal output p-norm 
is achieved with a family of input states of the form 
V T (8> V\ipe) with V unitary and \ipe) given by (|16p with 

» = «n- X (T^T=^)e(0,i)- 

ii) For p > fi c the optimal output 2-norm is achieved 
if and only if the input is the maximally entangled state 

IA>>. 

Moreover, the same conclusions hold for the minimal out- 
put entropy. 



In Section IIII CI this conjecture was proved for p = 
2. Since ||7||oo is simply the largest eigenvalue of 7, the 
observations at the start of Section IIV CI imply that the 
optimal input is achieved for a = 1 when p = 00. The 
conjecture then follows from the results in Section [Til Dl 

Theorem 2 Conjecture^ holds for p — 2 and for p = 00. 

Conjecture[T]is supported by extensive numerical work, 
as discussed in Section IIII El Additional evidence for the 
conjecture can be summarized as follows. 

a) When I < p < 2 and fi > \Z\i either Theorem Q] 
holds or the optimal input is achieved with a state 
of the form U <8> V\ipe) with U ^ V T and < 8 < 
sin -1 i. The latter would be unexpected, but has 
not been excluded. 

b) When p > 2 and fi < either Theorem [T] holds 
or the optimal input is achieved with a state of the 
form U <Z) V\ip e ) with U ^ V T and 6 > shr 1 \. 
Again, this seems unlikely. 

c) The heuristic argument described in Section HV Dl 
makes other behavior unlikely, especially for p > 2. 

d) For I < p < 2, Corollary 2] below excludes the 
possibility that the optimal input is maximally en- 
tangled for p < \i c . 

The conjectured behavior is very different from that 
for the shifted depolarizing channel on M4 which maps 

R^(l-fi)[XR+(l-X)(TiR)\h] +n(TrR)\0)(flflS) 

where \(3) is fixed and the quantity in square brackets [ | 
is easily recognized as the usual depolarizing channel on 
A/4. In this case the optimal input is always achieved 
using the vector \(3) which defines the shift. Changing 
from the usual depolarizing channel on M4 to a product 
of qubit channels ^a^^A dramatically changes the effect 
of the correlation introduced by a maximally entangled 
shift |/3) for values of /i < p c . 

Even without proving the conjecture, the next theorem 
implies that for I < p < 2 and for the minimal output 
entropy, the optimal input is never achieved with a maxi- 
mally entangled state when (i < p c . Moreover, to extend 
this result to p > 2, it would suffice to show that the 
optimal output is achieved with c v = 1. 

Theorem 3 For 1 < p < 2, the maximal output p-norm 
of the channel $/3 0!jLl ,A defined in JJ]) is achieved with a 
maximally entangled input if and only if that input is 
|/9o). The same result holds for the minimal output en- 
tropy. 

Proof: for I < p < 2 the optimal output is achieved 
with c v — 1. For any p the optimal input is maximally 
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entangled if and only if 6 = § in pT|). Then sin# = 1 
and (|31[) implies 

A-(l + A 2 )/ 

/ +M>| 2 M M V|a| 2 (l - M 2 ) \ 

-2A 2 

\M M Vl«l 2 (l - W 2 ) -2A 2 + M(l - |a| 2 )/ 

(44) 

from which it follows immediately that one eigenvalue of 
A is 1 + A 2 — 2A 2 = 1 — A 2 and the other two are 

1 + A 2 + ±M M ± ^ IM2 + 4A 4 - 2A 2 Af M (l - 2|a| 2 ) (45) 

which gives the largest p-norm when the term under the 
square root is largest, i.e., when |a| 2 = 1. But then, the 
results of Section [Til PI apply. QED 

Corollary 4 For 1 < p < 2 and \x < )i c the optimal 
output of the channel $/3 0iMi a is never attained with a 
maximally entangled input. The same result holds for 
the minimal output entropy. 

B. Majorization and trumping 

Let IV'opt) be the canonical vector for the family 
of input states which achieves the optimal 2-norm output 
for <&/3 0i A,/i> or, equivalently a state of the form U<g>V T |V>0) 
as described in part (i) of Conjecture [TJ If true, this 
conjecture would imply that for any p > 1, the output 
p-norms satisfy 

||$A,,A,M(l^opt)(?M)|| P > ||* j 9b,A,,.(l^>^l)llp(46) 

for any other input state One might, there- 

fore, expect that the eigenvalues of the matrix 
< i ) /3o,A, At (|i/ , opt)(V'opt|) majorize those of any other output. 
This is false, as can be seen from the following two ex- 
amples. 

First, let A = | and [i = ^. In this case /i > fi c = ^ 
so that the optimal input is \/3q), whose output eigenval- 
ues are {0,667,0.111,0.111,0.111}. A numerical search 
found that the input product state of the form (|A2D 
with ao = <22 = and 03 = = —ia\ yields an 
output state with eigenvalues {0.611,0.222,0.111,0.056} 
which is clearly not majorized by those of the preced- 
ing state. For an example with /i < /i c , consider A = 5 
and fi = j. In this case the optimal input gives an out- 
put with eigenvalues {0.596,0.141,0.141,0.123}. How- 
ever, the same product state now yields an output with 
eigenvalues {0.422,0.391,0.141,0.047} which, again, are 
not majorized by those of the preceding state. 

However, weak majorization x -< w y (which does not 
require equal 1 norms) is well-known to be equivalent to 
the stronger condition that < \\y\\ for all unitarily 
invariant norms (3ll . Corollary 3.5.9]. This is known as 
the Ky Fan dominance theorem [32]. Recently, Aubrun 



and Nechita [1^, [3(| proved an analogous result for £ p 
norms with p > 1 in which majorization is replaced by 
a relation known as "trumping" . This concept, intro- 
duced in connection with the phenomenon known as en- 
tanglement catalysis [28[ plays an important role in quan- 
tum information theory. If there is a vector z such that 
x (8> z -< y (& z, then one says that the vector y trumps x 
and writes a; <t V- Moreover, it follows that ||x|| p < \\y\\ p 
for all p > 1. Aubrun and Nechita [29|, [3(| showed that 
\\ x \\p < \\y\\p f° r ah P > 1 if and only if x is in the closure 
of the vectors trumped by y. Here, the closure is taken in 
the l\ norm and includes arbitrarily large catalyst vec- 
tors z. For a precise statement and other formulations, 
see 0,113. 

This leads us to the following reformulation of Conjec- 
ture [1] 

Conjecture 5 Let ^p ^.\ be a channel on M4 as defined 
in ([7]) and let fi c = |Et2 ■ 

i) For < fi < fi c , let ip opt be as in (j4l)|) . Then 
if W) 7^ ® ^OlV'opt} for some unitary V , the 
eigenvalues of $>p 0l \ tl j, (\ip)(i^\) yield a vector in the clo- 
sure of the set of vectors trumped by the eigenvalues of 

$/3o,A,/ i (|V ; opt)(^opt|)- 

ii) For (i > fi c and any input \^p) 7^ |/3o) the eigen- 
values of ^/3 0> A,At yield a vector in the closure 
of the set of vectors trumped by the eigenvalues of 
^ ,a,m(I/3o)</3o|). 

Conjecture [1] would follow immediately from Conjec- 
ture O It is, therefore, tempting to seek an analytic 
proof of this conjecture by seeking a catalyst z. Unfor- 
tunately, this is not easy in general. In our case, the fact 
that (1 — — A 2 ) is always an output eigenvalue, al- 
lowed a reduction to an effective 3-dimensional problem. 
For d = 3, it is known [28j that trumping with a finite 
dimensional catalyst z is never possible, although exam- 
ples are known [2a . [27j which use a infinite-dimensional 
catalyst. Thus, although we believe Conjecture O holds, 
it seems more likely to be established by proving Conjec- 
ture [T] and then applying Aubrun and Nechita's results 
[1^, [30] than by finding a catalyst. 

We can make additional reformulations suggested by 
the arguments in Section IIV1 For example, let A be 
given by (f3"Tj) with fi,\,0,<p fixed, and let v(|a|) denote its 
eigenvalues. Then if v(|a|) <t v (|1|) we could conclude 
that the optimal p-norm of ^p .^.\ is attained when |a| = 
1. The validity of Conjecture Q] would then follow from 
the arguments in Section ITlI Dl 
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APPENDIX A: ENTANGLEMENT 
PARAMETRIZATION 

Let \j3 ) = ^(|00) + |11}), and define = (/ ® 

<7/c)|/3o) for k — 1,2,3 and a k denotes the Pauli ma- 
trices. These four orthogonal maximally entangled 
states form an orthonormal basis for C 4 . The prop- 
erty Tri|/3fc)(/3 m | = \(JkVnn facilitates computation of 
reduced density matrices and entanglement of pure states 
represented in this basis. Our notation differs by a factor 
of i from the so-called "magic basis" introduced in [33| . 
When Ufc acts on the first qubit, 

= {ha> ll? < A1 > 

An arbitrary pure state € C4 can be written in the 
form 

3 

|V) = 5>*|ft) = (I®A)\0o) (A2) 

k=Q 

where A — ^2 k a k a k is in Mi- This gives an isomorphism 
between vectors in C4 and matrices in M^. We hence- 
forth restrict ourselves to normalized vectors for which 
IIV-II 2 = ELoM 2 = 1 = |TrAtA Combining 
with the observation that 02 1— > —a.2 takes A 1— > A T , we 
recover the well-known result that 

[iP) = (I ® A)\(3 ) = (A T ® I)\[3 Q } (A3) 

Since local unitary transformations do not affect the 
entanglement of the state in (|A2|) has the same en- 
tanglement as any state of the form 

(Z7®F)h/>) = 5^a fc (tf®J)(J®V)(J«>o*)|A>} 

k 

= (I®V)Y,ak{I®<r k ){I®U T )\{} Q ) 

k 

= (I®VAU T )\fj ). (A4) 

Moreover, by the singular value decomposition (SVD), 
one can choose U, V so that VAU T is diagonal and pos- 
itive, i.e., VAU T = D = aol + 030-3 > 0. By including 
a suitable permutation in U, V one can further require 
that the singular values arc in decreasing order, which is 
equivalent to ao > 03 > 0. Thus, we can choose U, V 
unitary so that 

|$) = (U®V)\i>) = (I®VAU T ) = a |/3o)+a« 



with ao > 03 > and aft = a| = 1. When this is rewrit- 
ten as in (jTTJ) , this implies that cos | > sin | > which 
impli es | e [0, j } and hence 9 G [Q, §]■ We emphasize 
that (jA5J) is not the most general pure state of this "diag- 
onal" form; instead V) should be viewed as a canonical 
state from which all others with the same entanglement 
can be written as (W ® V')\"ip). We use the same no- 
tation as in (fT3")l because the equivalence class of states 
defined by (fT2"]l corresponds exactly to the equivalence 
class of matrices with the same singular values. 

In the form (| A5|) it is easy to find the entanglement of 
\ip) because 

Tn|/3 )(/33| = Th|/9b)(A>|(/®<r3) = {\h)<rz = W A6 ) 

Therefore, the reduced density matrix of (|A5|) is 

7 = i(/ + 2aoa 3 CT3) , (A7) 

and its entanglement is 5(7) = /i(2aoa3) where /i(x) is 
the binary entropy 

h(x) = -i^Elnl±2 - i=3ini=2 . (A8) 

It this paper we consider subfamilies of states of the 
form CSS) with l/l/ T fixed. When UV T = I,U = V and 
the SVD becomes VAV> , which is the standard form for 
diagonalizing a self-adjoint matrix. Note that A is self- 
adjoint if and only if all a k in (|A2[) are real and this prop- 
erty is preserved by the transformation VAV^. Thus, if 
we begin with a state \ip) of the diagonal form (|A5|) with 
matrix D = aol + 03(73, then the family of states 

|^> = (V T ® 7*)$) = (/ ® WW^A,) = (I® A)|/3^A9) 

with V unitary all have real coefficients afc when written 
in the form ([X2j) . 

Since any maximally entangled state can be written as 
1/3) = (f/i ® C/ 2 )|/3o) = (/® U 2 Ul)\/3o), the vector |/3 ) 
can be regarded as the canonical representative of the 
equivalence class of maximally entangled states. This 
also gives a one-to-one correspondence between maxi- 
mally entangled states and unitary matrices in M2 via 
the relation 

|/3) = (I®U)\0o) = (U T <S,I)\(3o)- (A10) 

A matrix is unitary if and only if it can be written (up 
to an overall phase factor) as U = uqI + i J"), u k a k with 
Uk real. Thus, a state written in the form (|A2|) with ao 
real is maximally entangled if and only if Re a k = for 
fc = l,2,3. 

APPENDIX B: COMPARISON LEMMAS 

For n > Ti > r% it is elementary to verify that 

1 1 1 _ 

(n - r2)(ri - rs) (r-2 - n)(r2 - rs) (r3 - n)(r-3 - r'2) 
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and 



r 2 



(ri -r 2 )(r-i - r 3 ) (r 2 - n)(r 2 - r 3 ) (r 3 - ri)(r 3 — r 2 ) 





(B2) 



Lemma 6 Let P(x) — —{x — ri)(x — r 2 )(x — r 3 ) be the 
cubic polynomial with real roots n > r 2 > r 3 and let 
Q(x) = P(x) + S1X + S2 with 81,62 > 0. Then the roots of 
Q(x) are approximately s k w r k + ( rm - rk )( rn - rk ) ( r k8i + 
82) with k,m,n distinct. 

Proof: Near each of the roots, the tangent to P(x) is 

y = (x - r k )P'{r k ) = -{x - r k ) (t . m _ r J (r „_ rfc) (B3) 
from which it follows that near x = r k 

Q{x) » -(x - r k )(r rn - r k )(r n - r k ) + Six + S 2 . (B4) 
It then follows that Q(s k ) « if 

1 <5 2 



r k 



1 - 



r k + r k 



{r m -rk){r n —Tu)Si 

Si 



(r m - r k ){r„ - r k ) - 81 



(r m - r k )(r n - r k ) 
82 



(r m ~ r k ){r n - r k ) - 81 



QED 



(B5) 



Lemma 7 Let v be a vector in R3 with all 
vi > V2 > W3 > and let v k 1— > u>fc = 

+ (^-^ fc )V^-^ fc ) et,fc - Then H w ii x = i^H 1 anrf 

ll w llp > ll v llp w/ien e > and p > 1. 



Proof: The fact that ||w||i 
For p > 1 , observe that 



Ivlli follows from (IB2 



k k 



Thus, up to 0(e 2 ), 



H|g-||v||g 



qj 



- V2)(V1 - V3) (yi - V%){V2 - V3) 



(vi - u 3 )(u 2 — i> 3 ) 

1 



ep- 



Vl — V3 Wx — V 2 V 2 — Vz 



ep 



1 



-or 1 --?- 1 ) > 



> (B7) 
(B8) 

vi — V3 

where we used (|B1[) . and then the mean value theorem 
to obtain (|B8p with V\ > V\ > «2 and «2 > ^3 > V3. The 
inequality on the right in (|B7p follows from the fact that 



/(x) = x p is convex for p > 1. Although this suffices 
to prove the Lemma, the expression (|B8j) . will be useful 
when we need to compare competing terms. QED 
Lemma 8 Let v be a vector in R3 with all 
vi > V2 > v 3 > and let v k 1— » Wfe = 



e > 



-v k ){v n 



> 
< 



Then ||w||i = 

for p > 2 
for 1 < p < 2 



and /or 

(B9a) 
(B9b) 



Proof: The proof is identical to that of Lemma [7] above 
except that (|B7|) becomes 



1 



= 

V\ - v 3 

= ep(p- 1) 



mP- 1 



p-1 



i'2 

p-2 



W 2 - t>3 



,',P- 2 \ 



Vi - V 3 



(BIO) 



where, as before vi > v\ > v 2 and V2 > v 3 > v 3 . When 
p > 2, the function is convex; however, when 1 < 

p < 2, it is concave. Moreover, when p — 2 < 0, the 
expression on the right in (1B10[) is negative because v 3 < 
v 2 <vi- QED 







i(l - „)(1 + A 2 ) + § + 1^/^ + (1 - ^) 2 A 2 

*(1- M )(1-A») 

i(l- M )(l-A 2 ) 


|(1- M )(1 + 3A 2 ) + M 
-A 2 ) 

i(l- M )(l-A 2 ) 
1(1-^(1 -A 2 ) 


1(1 _ M)( i + a 2 ) + f - iVi^ + (1 - ^) 2A2 



TABLE I: Optimal output eigenvalues for 



APPENDIX C: EIGENVALUE COMPARISON 

The eigenvalues of A are easy to find for the two 
boundary cases of a = and a = 1. The lists in Ta- 
ble |TT] arc intended to be in decreasing order, but the two 
smallest eigenvalues may switch for very small M^. Al- 
though the largest eigenvalue is always greater for a = 1 
than for a — 0, the same is also true for the smallest. 
This precludes using majorization to conclude that all 
p-norms for a — 1 exceeds that for a = 0. However, as 
discussed in Section IV Bl the eigenvalues for a = 1 could 
still trump those for any a < 1. If so, the optimal output 
eigenvalues for p — 2 given in Table [J are conjectured to 
also be optimal for all p > 1 . 
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a = 






a = 1 




1 + A" + ^A 2 


-(1 


- A 2 ),? 2 / 


' 1 + A 2 + |M M + ^ 


j\Ml + 4A 2 - (1 - A 2 )S 2 + M M |AS*| 


1 - A 2 + M p 




\ 


1-A 2 




1 + A 2 - y^A 2 


-(1 


- A 2 )S 2 , 


' 1 + A 2 + \M„ - y 


/ \Ml + 4A 2 - (1 - A 2 )S 2 + Af M |AS| 



TABLE II: Eigenvalues of A with arrows showing expected 
increase and decrease with a. 
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